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We study cosmological perturbations generated from quantum fluctuations in multi-field infia- 
tionary scenarios in generalized Einstein theories, taking both adiabatic and isocurvature modes 
into account. In the slow-roll approximation, explicit closed-form long-wave solutions for field and 
metric perturbations are obtained by the analysis in the Einstein frame. Since the evolution of fluc- 
tuations depends on specific gravity theories, we make detailed investigations based on analytic and 
numerical approaches in four generalized Einstein theories: the Jordan-Brans-Dicke (JBD) theory, 
the Einstein gravity with a non-minimally coupled scalar field, the higher-dimensional Kaluza-Klein 
W), theory, and the R + R? theory with a non-minimally coupled scalar field. We find that solutions 

3 ' obtained in the slow-roll approximation show good agreement with full numerical results except 

■^C ' around the end of inflation. Due to the presence of isocurvature perturbations, the gravitational 

potential $ and the curvature perturbations TZ and C, do not remain constant on super-horizon 
scales. In particular, we find that negative non-minimal coupling can lead to strong enhancement of 
TZ in both the Einstein and higher derivative gravity, in which case it is difficult to unambiguously 
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I. INTRODUCTION 

^ , The beauty of the inflationary paradigm is that it both 1) explains why the present-day Universe is approximately 

\—^ ■ homogeneous, isotropic and spatially flat, so that it may be described by the Friedmann-Robertson- Walker (FRW) 
c/3 , model in the zero approximation, lit] and 2) makes detailed quantitative predictions about small deviations from 
Cd ' homogeneity and isotropy including density perturbations which produce gravitationally bound objects (such as 
^ galaxies, quasars, etc.) and the large-scale structure of the Universe. It is the latter predictions that make it possible 
to test and falsify this paradigm (for each its concrete realization) like any other scientific hypothesis. Fortunately, 
all existing and constantly accumulating data, instead of falsifying, continue to confirm these predictions (within 
observational errors). Historically, among models of inflation making use of a scalar field (called an inflaton), the 
original model or the first-order phase transition model B] failed due to the graceful exit problem, which was taken 
over by the new [pj and the chaotic B inflation scenarios where the inflaton is slowly rolling during the whole de 
Sitter (inflationary) stage. The latter property was shared by the alternative scenario with higher-derivative quantum 
gravity corrections pi (where the role of an inflaton is played by the Ricci scalar R) just from the beginning. Note 
that a simplifled version of this scenario - the R -\- E? model - was even shown to be mathematically equivalent to 
some specific version of the chaotic scenario Q (see also a review in Q). 

Turning to inhomogeneous perturbations on the FRW background, the inflationary paradigm generically predicts 
two kinds of them: scalar perturbations and tensor ones (gravitational waves) which are generated from quantum- 
gravitational fluctuations of the inflaton flcld and the gravitational flcld respectively during an inflationary quasi - de 
Sitter) stage in the early Universe. The spectrum of tensor perturbations generated during inflation was flrst derived 
in B , while the correct expression for the spectrum of scalar perturbations after the end of inflation was obtained in 
B. For completeness, one should also cite two papers pO| and pd| where two important intermediate steps on the 
way to the right flnal answer for scalar perturbations were made, in particular, in the latter paper the spectrum of 
scalar perturbations during inflation was calculated for the Starobinsky inflationary model |^. In order to obtain a 
small enough amplitude of density perturbations in all the above mentioned slow-roll inflationary models, the inflaton 
should be extremely weakly coupled to other fields. It is therefore not easy to find sound motivations to have such a 



scalar field in particle physics (see, however, 12 1^) 



Reflecting such a situation, the extended inflation [Q scenario was proposed to revive a GUT Higgs field as the 
inflaton by adopting non-Einstein gravity theories. Although the first version of the inflation model, which considers a 
first-order phase transition in the Jordan-Brans-Dicke ( JBD) theory ||l^ , resulted in failure again due to the graceful- 
exit problem |16|] , it triggered further study of more generic class of inflation models in non-Einstein theories |17| , in 
particular extended chaotic inflation |jl^ where both the inflaton and the Brans-Dicke scalar fields arc in the slow 
rolling regime during infiation. Note that the natural source of Brans-Dicke-like theories of gravity is the low-energy 
limit of the superstring theory ||l^,|22l with the Brans-Dicke scalar being the dilaton. 

Several analyses have been done on the density perturbations produced in extended new or chaotic infiation models 



|pi|-p5[, all of which made use of the constancy of the gauge- invariant quantity ( Pq], or its equivalent TZ [10 27| on 
super-horizon scales, and matched them directly to quantum field fiuctuations at the moment of horizon crossing which 
would be the correct procedure in a single component infiationary model. However, in the presence of two sources of 
quantum fiuctuations (i.e., the infiaton and the Brans-Dicke scalar field), C does not remain constant during infiation 
due to the appearance of isocurvature perturbations p8| . In such a case, mixing between adiabatic and isocurvature 
perturbations may occur due to ambiguity in the definition of the latter ones (see the discussion in the Sec. Ill A 
below) . 

Note that it is nothing unusual nor unexpected about non-conservation of the quantities ( and TZ even in the long- 
wave, or super-horizon, k <C aH limit since they are not conserved integrals of motion even in the one-field, k <^ aH 
case, if we follow the meaning of this term used in classical mechanics. Here, a{t) is the FRW scale factor, H = a/ a, 
where a dot denotes the time derivative, and k — |k| is the conserved covariant momentum of a perturbation Fourier 
mode. Namely, the "conservation" of ^ in the latter case is restricted to a part of initial conditions for perturbations 
for which the decaying mode is not strongly dominating. In the opposite case, since then C = 0(fc^$), the k^ term 



in equations for perturbations, e.g. in Eq. (3.16) below, may not be neglected even in the k <C aH case. Of course, 
any classical dynamical system with N degrees of freedom has exactly 2N conserved combinations of its coordinates 
and conjugate momenta irrespective of the fact if it is integrable or not (or even chaotic), but the functional form 
of these combinations is generically not universal and strongly depends on initial conditions. That is why no special 
attention is paid to such constant quantities in mechanics. The quantity C is just the examples of such a conserved 
combination for the growing mode. For the purely decaying mode, another conserved combination may be introduced 
(see Sec. HI below). It is the specifics of the inflationary scenario where the case of strongly dominated decaying mode 
is excluded that leads to the impression of the universal conservation of C in the one-field, k ^ aH case. Of course 
when more modes for each k appear in a multi-field case, generic non-conservation of ( becomes more transparent. 
This general remark explains numerous findings of non-conservation of this quantity in particular cases. However, 
this circumstance does not affect the predictive power of the infiationary paradigm at all since metric perturbations, 
in particular the gravitational potential $, may be calculated during and after infiation without any reference to ^ 
conservation or non-conservation. The only problem is that evolution of isocurvature modes of scalar perturbations (in 
contrast to adiabatic ones) is not universal, and its knowledge requires some additional assumptions about behavior of 
matter after the end of inflation (in particular, isocurvature modes disappear completely if the total thermodynamic 
equilibrium is reached at some moment of time) . 

A very important point in the derivation of spectra of primordial perturbations generated in the inflationary scenario 
is played by exact solutions of perturbations equations in the long-wave, or super-horizon limit A: — > 0. It is these 
solutions that give us a possibility to match quantum inflaton perturbations inside the de Sitter horizon during 
inflation to scalar perturbations during matter or radiation dominated eras, bypassing the study of physical processes 
in the Universe in between. That is why such solutions were sought and found for more and more complicated cases 



in numerous papers. In particular, in |29|, two of the present authors made the first correct analysis of the issue 
in the case of slow-roll inflation in the original Brans-Dicke gravity, extending the method used in [ pOpl] ] to find 
spectra of all modes of adiabatic and isocurvature fluctuations. This was further used to constrain the parameters of 
general scalar-tensor theories (e.g., the Damour-Nordvedt model p^) from the spectrum of the CMB anisotropy in 
p3| , p4| . General analytic formula for evaluating the spectral index in multi-field inflation was developed in |3^] , which 
neglected isocurvature modes, fc ^ solutions for all modes in the case of a factorizable potential U = Vi{Lpi)V2{(p2) 
were found in [ pq] . Finally, a general fc ^ solution for perturbations in the two-field case with an arbitrary potential 
U{(pi, (P2) was found in [|37| in the form of a functional over an infiationary background solution, from which explicit 
solutions for perturbations can be obtained in any case when the background solution may be explicitly integrated in 
the slow-roll approximation. 

However, there exists an old and completely different way to obtain super-horizon (fc -^ 0) solutions for perturba- 
tions, even without writing corresponding equations for them: the Lagrange method of variation of the background 
FRW solution a(i),0„(i),n = l,...iV with respect to all constants entering in it (here n numerates different scalar 
fields). In the case of matter in the form of classical fields, this method always directly produces all 2N physical 
modes of the k —^ solution for some quantities (in particular, for some gauge-invariant quantities, too), though to 
obtain the full fc — > solution for all gauge-invariant quantities one has additionally either to use the — i Einstein 



equations, or to integrate the i — j {i ^ j),i,j = li2,3 Einstein equations.^] The latter necessity is related to the 
fact that all gauge-invariant quantities constructed from a space-time metric and its first derivatives (including a 
Newtonian potential $) are necessarily non-local, in accordance with the Einstein equivalence principle. As a result, 
their expressions in terms of quantities which may be defined and measured locally contain k^ in the denominator. 
So, the next order of the series in powers of fc^ (beyond that which follows from the Lagrange method) has to be 
computed for their determination. 

Since the background FRW space-time is in the synchronous form, the Lagrangian method yields fc — ^ per- 
turbations in the synchronous gauge, too. Of course, this method produces solutions in the explicit from only if 
the background solution is known explicitly, too. Thus, its power in the two-field case is the same as that of the 
Mukhanov-Steinhardt functional. On the other hand, it is applicable for any number of scalar fields and during any 
stage of the Universe evolution, not necessarily at slow-roll inflation. 

In the case of adiabatic modes of perturbations (both growing and decaying ones), this method was not only known 
for many years, but has been already used by Lifshitz and Khalatnikov in a more advanced form - to produce non- 
homogeneous solutions near singularity without assuming inhomogeneous perturbations to be small. This is achieved 
by taking integration constants in the background FRW solution as arbitrary functions of spatial coordinates. That 
was how their quasi-isotropic solution |pq|, or the so called 7- functional solution 139] were constructed. In the case of 
isocurvature modes of perturbations, this method was recently considered in details in [[40|-[42|. 

We will use this method in the case of adiabatic modes of perturbations, since it is especially simple in this case. 
Also, it provides a simple reason for the universal constancy of the properly deflned growing adiabatic mode at super- 
horizon scales (different from that recently proposed in Q) which does not depend on any local physical process, in 
particular, on presence or absence of preheating. On the other hand, since the Lagrange method and the Mukhanov- 
Steinhardt functional are equally powerful in the two-field case in the slow-roll approximation (so that we may use 
any of them), we will use the latter formula to derive explicit expressions for non-decreasing isocurvature modes of 
perturbations. 

At present there are many generalized Einstein theories which can provide infiationary solutions, e.g, generalized 
scalar-tensor theory, Einstein gravity with a non-minimally coupled scalar field, higher-dimensional Kaluza-Klein 
theories, f{R) gravity theories. Making use of the conformal equivalence between these theories, cosmological pertur- 
bations can be analyzed in a unified manner (see, e.g., [|44| ). In the present paper we analyze density perturbations 
generated in slow-roll inflationary models for a general class of generalized Einstein theories in the presence of two 
scalar fields. We make use of the conformal transformation |}45|| which transforms the original, or the Jordan frame 
to the Einstein frame in which equati ons are somewhat simpler. In particular, the i — j {i ^ j) Einstein equations 
directly lead to the universal relation ( |3.§| ) for all models considered in our paper. If this class of gravity theories is 
considered as a low-energy limit of the supcrstring theory, then the Jordan frame is also called the string frame. 

Though it is always possible to write fc — > solutions for perturbations in a functional form, actual evolution of 
perturbations depends on specific gravity theories. In the JBD theory, where the Brans-Dicke parameter is constrained 
to be w > 3500 from observations Q , the gravitational potential is dominated by adiabatic perturbations |2^j3J] , in 
which case the variation of TZ is restricted to be small as we will see later. On the other hand, it was recently found 
that negative non-minimal coupling with a second scalar field other than infiaton can lead to significant growth of 
C by the analysis in Jordan frame [}47| . In this case it is not obvious whether the slow-roll analysis provides correct 
amplitudes of field and metric perturbations, since these exhibit strong enhancement during inflation by negative 
instability. In this paper we will investigate the validity of slow-roll approximations by numerical simulations in the 
Einstein frame. We will also work on the evolution of cosmological perturbations in a higher-dimensional theory and 
the R + B? theory with a non-minimally coupled scalar field. The latter corresponds to the case where explicit and 
closed forms of solutions in the super- horizon limit are obtained by slow-roll analysis in spite of a coupled form of the 
effective potential. 

The rest of the present article is organized as follows. In Sec. II we present the Lagrangian in the Einstein frame and 
introduce several generalized Einstein theories which can be recasted to the Lagrangian by conformal transformations. 
Then in Sec. Ill basic equations and closed form solutions for super-horizon perturbations are given. In Sec. IV we 
apply the results of Sec. Ill to specific gravity theories, namely, the JBD theory, the Einstein gravity with a non- 
minimally coupled scalar field, the higher-dimensional Kaluza-Klein theory and the i?^ + (l/2)^i?x^ theory. In order to 
confirm analytic estimates, we also show numerical results by solving full equations of motion. We present conclusions 
and discussions in the final section. 



'^For matter in the form of A'^ hydrodynamic fluids, only A'^ non-decreasing modes, including the most important growing 
adiabatic mode, may be obtained in this way. 



II. INFLATION IN GENERALIZED EINSTEIN THEORIES 



Consider the following two-field model with scalar fields ipi and (p2'- 



(2.1) 



where k^/Stt = G is the Newton's gravitational constant, F(ipi) is a function of (pi, and U{(pi,(p2) is a potential of 
scalar fields. Many of the generalized Einstein theories are reduced to the Lagrangian ( p.l[ ) via conformal transfor- 
mations |45|. We have the following theories, which may provide inflationary solutions. 



1 . Theories with a scalar field ^ coupled to gravity whose action is written by 



S ~ I d^x\J —g 



fWR{g) - h{i;)(y^)^ - i(V(/.)2 - V{^) 



(2.2) 



where V{(l>) is a potential of inflaton, 0. In this work we consider the following theories. 
(a) Jordan-Brans-Dicke (JBD) theory with a Brans-Dicke field, ip [|5|. In this case / and h are 

f = — /i= ^ 
167r' 16mp' 



(2.3) 



where ui is the Brans-Dicke parameter which is restricted as lj > 3500 from observations pq]. Making a 
conformal transformation, 



QfJ-V — ^' 9tJ.V: 



where 



n' 



k2 / HX 

— w = exp — . 

Stt \V^ + 3/2 



we obtain the action in the Einstein frame (2.1) with replacement, 

'Pi ^ X, V2^ (/), 
and 

F=(/3/4)^X, U{x,cj)) = e-''-W{<f,), 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



with P = ^87(2^7+3). 
(b) Non- minimally coupled massless scalar field, ^, with an interaction, {\/2)S,RijJ^ [EtIEs]. In this case / and 



h read 



/ 



1 - C/^V^ 
2k2 



^=2- 



(2.8) 



Applying the conformal transformation ( ^.4| ) with VP = 1 — ^n'^ip'^ and defining a new field x in order for 
the kinetic term to be canonical as 



'l- (1-60^^2^2 
(1-$k2^2)2 



dtp, 



we obtain the action (|2.1| ) with replacement (^_^) and 



^ = ^in|i-^-Vl, f^(x.^) = ^-'"^^^m-(T^^i|^- 



(2.9) 



(2.10) 



The induced gravity theory [^ is also described by the action, ( p.2[ ). In this theory the scalar field ip has its 
own potential of the form, ¥{1^) = (A/8)(V'^ - fff with / = (e/2)v/and h = 1/2. 



2. The higher-dimensional theories where the inflaton, 0, is introduced in N = D + A dimensions 

R 1 



S= I d"" 



xy/-g 



2k2 



-(V0)2-y(<^) 



(2.11) 



where /t^ and R are the A'^-dimensional gravitational constant and a scalar curvature, respectively. We compactify 
the A'^-dimensional spacetime into the four-dimensional spacetime and the _D-dimensional internal space with 
length scale, 6. Then the metric can be expressed as 



dsjq — cjfiydx'^dx^ + b dsjj, 



(2.12) 



where g^j/ is a four-dimensional metric. Assuming that extra dimensions are compactified on a torus which has 
zero curvature, |j one gets the following action after dimensional reduction |2l[ : 



S = / d'^Xy 



6\^ 1 



'^\boJ 2^2 



R + d{d 



_ 9^ ,ri 



62 



-(V<^)2-F(0) 



(2.13) 



where bo is the present value of 6, and R is the scalar curvature with respect to 5^,^. In order to obtain the 
Einstein-Hilbert action, we make the conformal transformation (p.4[) with a conformal factor, 



n^ = exp Z? 



where a new scalar field, x, is defined by 



X_ 
Xo 



X = Xoln 1 ^ ) ' '^ith Xo = 



(2.14) 



2k2 



1/2 



(2.15) 



Then the four-dimensional action in the Einstein frame can be described as (2.1) with replacement (2.6) and 
F = 0, C/(x,(/.) = exp(-/3/cx)F(0), with /3 = 



2D 

D + 2' 



(2.16) 



Note that when inflaton is introduced in the fo ur-dimension al action after compactification, we find F — 
g-(/3/2)«x and [/(x,^) = e-'^'^^V{(j)) with /3 = ^8D/{D + 2). In this case, however, we do not have infla- 
tionary solutions since the effective potential does not satisfy the condition: (3 < y/2, which is required for 
power-law inflation to occur (see the next section). 

3. The f{R) theories where the Lagrangian includes the higher-order curvature terms, i.e., df /dR depends on the 
scalar curvature R [H: 



S^ I d^ 



XV -5 



f{R)-l{^xr 



(2.17) 



In this case the conformal factor 



n^ = 2k2 



df 



dR 



(2.18) 



Note that there exist other methods of compactifications. One of them is the compactification on the sphere ||5C|| , in which 
case stability of extra dimensions and the evolution of cosmological perturbations during inflation are studied in [Rl|J52[]. 



describes a dynamical freedom in the Einstein-Hilbert action. Introducing a new scalar field 



7 In 



2k^ 



df 



dR 



(2.19) 



the action in the Einstein frame is described as ( p.l| ) with replacement, 



(2.20) 



and 



-F = — =, C/(0, x) = (sign) exp ( —kcJ) 



^i?(0,x)exp(f J-mx) 



(2.21) 



where sign — {df/dR)/\df/dR\. For example, in the R^ theory [g[ with a non-minimally coupled massless x 
field, i.e.. 



1 ^ ^2 1 



f(R) = -—R + qR'--^Rx\ 



(2.22) 



the effective two-field potential is described as 

,4 



^^^' ^) = J^'^^'^"'^^' (^^^^'^'^' - ' + ^"'^') ' 



(2.23) 



where we have chosen a positive sign. In this case the (j) field behaves as an inflaton and leads to an inflationary 
expansion of the Universe. 



III. COSMOLOGICAL PERTURBATIONS IN TWO-FIELD INFLATION 



A. Basic equations and solution for the adiabatic mode 



Let us first derive the background equations. Variation of the action ( |2.lD yields the following background equations 
for the cosmic expansion rate H — a/ a and homogeneous parts of scalar fields: 



1 



1 



i^^ = y^2^^ + --^^ + t/ 



(3.1) 



H 



Wi 



-2F ■2\ 

e V2) 



(3.2) 



(fi + ZHlpi + C/, 



VI 



F^.e-^'^ifl^O, 



(3.3) 



„2F 



^2 + 3Hif2 + e^-^ t/^^, - 2F^,<fiifi2 = 0, 



(3.4) 



where a prime denotes a derivative with respect to ipi. A generic solution of this system contains 5 arbitrary inte gration 
constants (two constants appear from the solution of Eq. (3.3), two - from Eq. (3.4), and one from Eq. (3.1), while 
Eq. (|3.2| ) is a consequence of the other equations). However, one of these constants corresponds to a trivial shift of 
the cosmic time t. The Lagrange variation with respect to these constant yields a gauge mode. So, variation with 
respect to only 4 constants may be used to produce physical solutions for perturbations in the long- wave limit. 

Moving to perturbations now, we restrict ourselves to the spatially flat FRW background, first, for simplicity and, 
second, because recent data on angular fluctuations of the cosmic microwave background (CMB) convincingly confirm 
the absence of any significant spatial curvature of the Universe within a few percent accuracy (see, e.g., |Q). Then 
a perturbed space-time metric has the following form for scalar perturbation in an arbitrary gauge: 



ds'^ = -{I + 2A)dt^ + 2a{t)B^,dx'dt + a'^{t)[{l + 2D)S,j +2E^,^j]dx'dx^ , ij = 1,2,3, (3.5) 

where a comma means usual flat space coordinate derivative and A is the flat 31? Laplacian. In the synchronous 
gauge, A ^ B = 0, D = {X + /i)/6, AE = —A/2 in the Lifshitz notations. In the longitudinal gauge, B — E — 
0, A — f^, D = — ^I/, and $,^ are gauge-invariant potentials [p4|-p6[. Further, we assume the exp(ikx) dependence 
for each Fourier mode k and omit the subscript k = |k| in expressions for time-dependent parts of perturbations. 
Note the useful relations between A, fj, and $, ^, and also between scalar field perturbations in the synchronous gauge 
5'fs,n, n = 1, 2 and in the longitudinal gauge Sipn (the latter quantities are gauge-invariant actually): 



* = -2fc2rf^K^)' '^ = —Q^^ + ^'^ + W^^ S^„=S^s,n-^X. (3.6) 

Another useful gauge-invariant scalar field perturbation is given by the Mukhanov-Sasaki variable ||57[|: 



g„ = 6ip„ + -^* = 6(ps.n - 777 (^ + /^) • (3-7) 

For all models considered in our paper, we have the following relation in the Einstein frame 

$ = *, (3.8) 

which follows from the i — j {i ^ j) Einstein equations taking into account the fact that anisotropic stresses vanish 
at the linear order there. Transforming back to the Jordan frame, this relation does not hold generically pqj44[ |. 

Now we can derive one solution for super-horizon perturbations using the Lagrange method without writing equa- 
tions for perturbations. This is possible since one of integration constants of a background FRW solution, namely. 



that which appears by integrating Eq. (3.1) is trivial: it is simply a multiplier oq of a(t). Of course, due to invariance 
of measurable quantities with respect to equal rescaling of all 3 spatial coordinates x* , this constant does not appear 
in variables like H{t) and (pn{t). So, in the fc -^ limit, 

/x = 6 =3h — const; X,dips n — Oik^h); q„ — ^h—^. (3.9) 

Go ' 2H 

This formula is valid both in the Jordan and the Einstein frames. By definition, this partial solution will be called 
the growing adiabatic mode (we will discuss the ambiguity of this definition later). It is clear that this solution exists 
for any form of the gravity Lagrangian and the matter energy-momentum tensor. In particular, presence of fast 
oscillations in a background solution does not affect it, too. The only things which are needed for its existence are 
the spatial flatness and isotropy of the background metric. That is why the derivation of the spectrum of adiabatic 
perturbations generated during inflation in the second of Refs. y, where this method was used even in a more general 
form (valid when \h\ is not necessarily small), does not depend on any physics between inflation and the present era. 
Note also that the invariance of measurable quantities with respect to different rescaling of spatial coordinates with 
the total volume fixed leads to the constancy of the non-decreasing (quasi-isotropic) mode of gravitational waves in 
the super-horizon regime. That property played a crucial role in the derivation of the spectrum of gravitational waves 
produced during inflation in B. 



Let us now calculate the gauge-invariant quantities $ and Sifn for the solution (3.9) in the Einstein frame where 
the relation (^^) holds. In the synchronous gauge, the latter relation reads 



2 



A + 3iJA--^(A + /i) =0 . (3.10) 



Thus, in the fc — > limit, we get 



X^^fadt, $ = * = -J (l-— / adt) , Sipn^~h^ I adt (3.11) 



"^ Jti ^ V u jti / 2a jt^ 

(ti depends on k generically). A shift in the ti produces one more super horizon solution - the decaying adiabatic 
mode: 

A = /lla"^ fi^O{hik'^), <^ = -^ = hiH/2a, S(pn = -hi(pn/2a . (3.12) 



Thus , the growing adiabatic mode ( |3.9D is defined up to an addition of some amount of the decaying adiabatic mode 
( |3.12 ). In the inflationary scenario, ii is the moment of the first Hubble radius crossing during inflation for each 



k, so there is n o ambiguity at all. Thus, as a whole the adiabatic solution (mode) is defined unambiguously by 



Eqs. ( |3.9| , 3.11 ). On the other hand, one may always add some amount of the adiabatic mode to other, isocurvature 
solutions (modes). So, the definition of isocurvature modes is not unique. In particular, this ambiguity may be used 
to make $ = ^I^ = for an isocurvature mode in the Einstein frame at some chosen moment of time, e.g., at the end 
of inflation or at the moment when the full thermal equilibrium is reached (if the latter occurs at all). Of course, this 
choice does not affect any observable quantities. 

A useful gauge-invariant quantity is the comoving curvature perturbation TZ |IQ,E^ : 

7^ = -i(A + ^) + -^(A + /i) = * - :^(* + i/$) . (3.13) 

A similar quantity is the curvature perturbation on uniform density hypersurfaces introduced in P6| : 

(^^D~H^-^=(-- -^] {\ + ^^)-^(, = -n+ -^* (3.14) 



where p is the total energy density of matter and 5p - its perturbation. For the solution (|3.9D , we get 

7^ = - C = -h/2 ^ const ^ 0. (3.15) 

Thus, we have the theorem: 

In the superhorizon limit k ^ Q and for H ^ identically, there always exists one solution for perturbations (the 

growing adiabatic mode) for which TZ = —C, is conserved in the leading order in k^ (apart from the vicinity of points 

where H — 0). 

Let us emphasize that this statement is valid both in the Jordan and the Einstein frames. 

Pathological behavior of TZ and C at the moments of time when H = where they diverge (if terms of the next 
order in k^ are taken into account) is solely an artifact of their definition. The potentials 4> and '5 remain regular and 
small near these po ints , so n o consideration of non-linear eff ects i s required there. Moreover, it immediately follows 
from the formulas (^^, |3.11 ) that the same constant value ( 3.15| ) of TZ and C is restored after passing through any 



point where H — 0. This behavior of ^ is clearly seen in our numerical calcul ation s in Figs. 1 and 3 below. 



What about conservation of TZ and C. for the decaying adiabatic mode ( 3.12 ) ? Here, there is a subtle point. 
In the leading order, TZ — ( ~ for this mode. However, one may not say that these quantities are conserved 
(even approximately). Real (approximate) conservation would require H\TZ/TZ\ <C 1 that is not valid if the decaying 
adiabatic mode is strongly dominating p8|. Also, the relation TZ « — C is not correct for this mode. On the other 
hand, in the Einstein frame one can introduce the gauge invariant quantity T = a^/H which does conserve in the 
super- horizon limit in the leading order in k'^. 

Finally, for other solutions (isocurvature modes) TZ and Care not conserved in the super-horizon limit, too. This 
can be easily seen by considering the time derivative of TZ ||3^,^ : 

where 



TZ=^ — <P + H[^-^]Z, (3.16) 



i^l + e'^F^ly 






The quantity 



Stp^ip^ = H I —, — ] , (3.18) 



represents a generalized entropy perturbation between ipi and (/?2 fields |59,M. In the multi-field case, Sip^ip^ ^ and 
Z ^ generically. So, the presence of isocurvature perturbations leads to the variation of 7?. ~ — C- 

Note also that we use the terms "adiabatic perturbations" and "isocurvature perturbations", as is commonly done, 
to denote different solutions (modes) of the same physical variables. This should be contrasted with the approach of 
the recent paper [p9| where "adiabatic" and "entropy" fields are introduced which are certain linear combinations of 
initial fields (pn{t), and then adi abatic and entropy perturbations mean perturbations of these fields. However, for 
our adiabatic solution (3^, |3.1l|) the perturbed entropy field 6s ex {(piSip2 — (p2S(pi) = in the super-horizon limit. 



even if the background field trajectory in the scalar field space is curved. Thus, the adiabatic mode is not sourced by 
isocurvature (entropy) modes, in contrast to results of |59f| . 

That is all what may be obtained without solving equations for perturbations. So, to proceed further, equations 
for time-dependent parts of metric and field fluctuations have to be written. They read: 



$ + 4iJ$ + K^U<^ = — [ipidipi - {U^^, + F^^,e-^^ip2')5ipi + e-^^ip2Sip2 - U,^Jip2] , (3.19) 



^ - if j * = -y [93i<5<^i + [iH^i + t/,^, - F^^.e^'^^ ^1)5^1 + e~^^ ^28^2 + [U,^, + ?,H^2e~^^)8^2\ , (3.20) 






(3.21) 



r/,2 
'2 



5Cpi + 3HSipi 



^^ + C/,yivi - (e ),<^i^iY 



Sipi + 2F^^^e '^^ip25^2 + U^^^^^Sip2 



(3.22) 



Sip2 + (3H - 2F)Sip2 +[^ + e'-^t/,^,^, ) Sip2 - 2FSipi + e^^ (2F^,t/^, + U^,^,) Sipi 
= 4(p2$-2e2^C/,^,$. 



(3.23) 



The relation (3.20) clearly indicates that metric perturbations are determined when the evolution of scalar fields are 
known. 



B. Closed form solutions in slow-roll approximations 



The use of the slow-roll approximation allows us to obtain closed form solutions for isocurvature perturbations in 
the long- wave limit p9 33 3qj37[|. Under this approximation, the background equations are simplified as 



H' = -^U, 



(3.24) 



3H(fi + C/,y, = 0, 



(3.25) 



SHip2 + e^^U,^, = 0. 



Combining Eqs. (^^)-( ]3.26| ) with Eq. (|1|), we find 



fi = 



k2 U 



V2 



k2 ij 



(3.26) 



(3.27) 



H 



1 

2^ 



U 



VI 



U 



2F I ^.^2 
U 



(3.28) 



In JBD and higher-dimensional theories where the potentials take the form, U{(pi, 1P2) — e '^'^^^V{(p2), it is straight- 
forward to show that the scale factor evolves as power-law. In this case integrating ipi = f3H/K over t, we find 



(pi[t) = -In h ipif = z + ipif, 



(3.29) 



where a subscript / denotes the value of each quantity at the end of inflation and z is the number of e-folds of 
inflationary expansion after the time t. Assuming that V{^<^2) takes a constant value Vo during inflation, one finds 



a = ao 



1 ^/^ B^ ' "" 



(3.30) 



where ao and to are constants. Then we have a power-law inflationary solution when /? < \/2. For example, in 
the JBD case with potential, (2.7), inflation takes place with a large power exponent because (3 is constrained to 



be /3 = yj'&l'luj + 3 < 0.034 from observations p3|. In higher-dimensional theories with potential (2.16), inflation is 

realized for arbitrary extra dimensions £>, because the condition, /? < -\/2, always holds. 

Let us consider large scale perturbations with fc ^ all . Neglecting $ and those terms which include second order 
time derivatives in Eqs. (3.21 )-(B^), one finds 



^=2H ('^I'^'^i 



-IF 



(^2(5(^2 



(3.31) 



3HSipi + U^^-^^^Sipi + U^ip^^^5ip2 + 2[/,^i$ = 0, 



(3.32) 



2,H5^2 + (e'^t/,^,)^^^ (5(^1 + (e^^f/,^,)^^^ 6^2 + 2e2^t/,^,$ = 



(3.33) 



Note that this approximation may not be always valid especially when perturbations exhibit nonadiabatic growth 
during infiation. We will check its validity in the next section. 

In the slow-roll approximation, the expression (3.11) for the growing adiabatic mode is simplified: 



$ = ^ = hJL , 5^^ ^ ^h^ 



(3.34) 



To fi nd a generi c sol ution, we introduce new variables, x and y, with Scpi = U^^-^x and Sip2 = e U^^^U- Then 



Eqs. ( |3.32D and (|3.33| ) yield 



U U e^^ 



(3.35) 



3Hy 






(x - y) + 2$ = 0. 



Subtracting Eq. (3.35) from Eq. ( ^.36 ), we obtain the following integrated solution 

A 



y - x == (^3 exp 



3H 



-dt 



where Qs is a constant, and 



Taking notice of the relation 






-2F 



u^. 



u 



V2 



(3.36) 



(3.37) 



(3.38) 



$ = 



2H 



Ux + t/,<p2<^2(y - x) 



K 

2H 



Uy + U^^^ipi{x-y) 



and making use of Eq. ( 3.37 ) and background equations ( ^.24 )-( 3.26| ), we find 



V3 

' U 

03 

' u 



TT TT TT p2F 

h (V,^2i^2 



U, 



J , 

—dt, 



H (e^^ t/, 



,C/.o 



y2Ayi ^,yi , tj ■ 



U 



J , 

Jjdt, 



(3.39) 

(3.40) 
(3.41) 
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where 



J =U cxp 



_A_ 
3H^ 



dt 



Then the final closed-form solutions for long-wave perturbations are expressed as 



&Vi = (lnC/),yi 



Q1 + Q3 (ln(lnC/),^ J ^ Jd(^: 



(3.42) 



(3.43) 



5ip2 = e-"" (InC/), 



1P2 



Q2-Q3J (ln(e2^(lnC/),^,))^yd(pi 



(3.44) 



with 



* = -2 [(lnC/),vi'5(^i + i\nU),^J^2] , 



J = exp 



t ^ 



(In (e2^(lnC/),^,))_^^ d(^i -H (ln(lnt/),^J^^^ dy.2 



(3.45) 



(3.46) 



Here integration constants, Qi, Q2, and Q3 satisfy the relation Q2 = Qi+ Qs, which comes from Eq. ( ^3.37 ). The Q3 
terms appear due to the presence of isocurvature perturbations. These constants are evaluated by the amplitudes of 
quantum fluctuations of scalar fields at horizon crossing, i*. The fluctuations are generated by small scale perturbations 
(fc > aH), so that they can be considered as free massless scalar fields which are described by independent random 
variables P9|j3l| ]. Then the field perturbations when they crossed the Hubble radius (fc ~ aH) are written in the form: 



^¥'i(k)|t 



HjU) 
^/2fc3 



evi(k), S(p2{'k)\t=f,, 



^^^^^(*-)e^,(k). 



/2fc3 



Here e^p^ and e^p^ are classical stochastic Gaussian quantities, described by 

(e^, (k)) - (e^, (k)) = 0, (e^, (k)e;^ (k')) = 6,,S^'^ (k - k') 



where i,j — ipi,ip2- From Eqs. ( |3.27| ), ( p.43D , and ( p.44D , we find 






K 



Q1 + Q3 (ln(lnC/),^J ^ Jd^2 



t. 



(3.47) 



(3.48) 



(3.49) 



5^2 

<P2 



K 
ll 



2 r 



1 + Q3-Q3 (ln(e2^(lnf/),^,))_^yd(pi 
Making use of Eqs. (3.47), (3.49|), and (3.50), the integration constants are expressed as 



Qi = - 



H^U) /e^,(k) 



^^^ H^t^) / e^,(k) ^pe^^jk)- 

K2^2fc3 \ (pi V>2 / t. 



K"V2fc3 V Vl / t. 

In the next section, we will investigate the evolution of large-scale perturbations in specific gravity theories. 



(3.50) 



(3.51) 



IV. APPLICATIONS TO SPECIFIC GRAVITY THEORIES 



Among the generalized Einstein theories which we presented in Sec. H, most of them take the following separated 
form of potentials except for the /(i?) theories: 



U{(pi,ip2) = Vi{ipi)V2{f2)- 

In this case, we have an integration constant by making use of Eq. (^.27|): 



(4.1) 
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C = -k" 



Vi 



e dipi + K^ 



^2 



1^5 



2,V2 



-dip2 



(4.2) 



This characterizes the trajectory in field space in two-field inflation. 

Since J = e^^^^"-*^*^ in the separated potential, (4.1), Eqs. ( 3.43| ) and ( 3.44| ) are easily integrated to give 



Vi V2 



together with the gravitational potential. 



^^-Uf)^.4(|)^Q.-.g3e-: 



(4.3) 



(4.4) 



Introducing new integration constants, Ci = —k^{Qi +Q3,e^^*) and C3 = —n^Qze^^' , and making use of Eq. ( 3.28 ) 
one finds 



,5^i--(Ci-C3) 






5^2--[Cie2^-C3(e^^-l)] 






(4.5) 



H 



where e^p-^ and e^^ ^-re given by 



1 fV{ 



2k^ V^i 



1 fV^ 



2k^ VV^2 



From Eq. ( 3.51 ) Ci and C3 are expressed as 



Ci 



V2fc3 



2F\ '--'P 



(1-e-) 



e<pi(k) 3_peip2(kj 



'/'I 



",52 



C3 



g^(^.) 
%/2fc3 



3-F£v2_W. _ 2F ^yi(^) 



V2 



•fl 



(4.6) 



(4.7) 



(4.8) 



The gravitational potential can also be decomposed in a different way. For example, let us introduce the integration 
constants Ci and 6*3, defined by 



Ci = -k'Qi 



^2g2(F.-F,) JJ2^^ 



l + a/ 



Q3 



/2P 



1- 



e2(^-^^)\ e^,(k) , e3F-2F,g^^(k) 



1 + a/ / (^1 1 + a/ <y52 



r^ — ^2„2(F,-Ff)^ -^ V-* 



2/c3 



^3F-2F 



'j 6^2 (^^) _ 2(F-F/) Syi(J^) 



¥52 



'/'I 



where the subscript / denotes the value at the end of inflation, and 



^V2 



Then Eq. (4.4) reads 



$ = -Ci-^-C3 



'^yi 



1 + Uf \ 1 + a/ 



=2F 



.2Ff 



(4.9) 



(4.10) 



(4.11) 



This decomposition corresponds to the case where the second term in the rhs of Eq. ( 4.11 ) vanishes at the end of 

inflation. 

During slow-roll we have that \H/H'^\ < 1 in Eq. ( |3.28D , which yields from Eqs. ( |3.13D , (jij), and ( |4.1l| ) 



as 



7^ 



H 



H 



-'^ = Ci-Cz'^^ 



e^, + (e^^ - l)e 



^Vl 



^V2 






(4.12) 
(4.13) 
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Note that both decomp ositi ons coincide each other in the hmit a/ — > and Ff -^ 0. When a/ and Ff are nonvanishing, 
the second term in Eq. (4^) gives contributions to the gravitational potential, $. When this term is negligible relative 
to —CiH/H^ during the whole stage of inflation, the first and second terms in Eq. (4^) can be identified as adiabatic 
and isocurvature modes, respectively. However, in some generalized Ein stein theories which we discuss in the following 
subsections, the final $ is dominated by the second term in Eq. (4.6). In those cases we can no longer regard the 
second term as the isocurvature mode at the end of inflation. 



The relation ( 4.12 ) or ( 4.13 ) indicates that TZ is not generally conserved during inflation when both fields are 
evolving due to the presence of isocurvature perturbations. The generalized entropy perturbations are written as 



S^,^,^-C,e-^^ = -C,^'-^'~^\ 



(4.14) 



Since S^p^^p^ and Z do not vanish generally, these work as a source term for the change of TZ. The evolution of 7?. 
depends on specific gravity theories as we will show below. 



A. JBD theory 



Let us first apply the results in the previous section to the JBD theory with Eq. 
( [4.12[) are 



7h. In this case Eqs. (4.6) and 



$ 



-Ci 






C-3 



r/?2 



+ (g(/3/2)«x _ 1)£^ 



(4.15) 



7^ = Ci - C3 



1- 



g(/3/2)«x + /32/(2e^) 



(4.16) 



In the JBD theory, /3 is required to be /3 < 0.034 from observational constraints, which yields e^ = /3^/2 ^ 1. 
In addition to this, the value of x should be practically vanishing, Xf — Oj ^t the end of inflation in order to 
reproduce the present value of the gravitational constant, because x generally grows only as logarithm of t after 
inflation and is even constant during radiation-dominated stage in this t heory . These lead to Ff = {f3/A)KXf — ^-^d 
af — e~'^^f exf/f-4>f — 0, which implies that the deco mposition of Eq. ( 4.11 ) looks almost the same form as that of 
(4.6). Since the second term in the rhs of Eq. ( 4.15| ) is negligible relative to the flrst term after inflation, the term 
proportional to Ci represents the growing adiabatic mode [p5| , while the one proportional to C3 corresponds to the 
isocurvature mode 129,641 • During inflation, TZ evolves due to the change of the term. 



^^g(/3/2)KX+^2/(2e^)_ 



(4.17) 



As a specific model of inflation, let us first consider chaotic infiation driven by a potential, V{4)) = A2„(?!)^"/(2ri). 
In this model, inflati on o ccurs a t la rge and it is terminated when \4>/ 4>\ becomes as large as H a.t (j) = (j)f ~ y/2n/K. 
Making use of Eqs. ( |3.29| ) and (U) we find 



V{ip) 






e-^^/^) ^ z. 



(4.18) 



The error in the last expression is less than 1.7% for z < 60 and /3 < 0.034. Then Eq. ( 4.17 ) is rewritten as 

W=(\--\e^P/^)-^ + -+'^. (4.19) 

\ nj n 2n 

Since W is constant for n = 2, TZ is conserved in the quartic potential, V{(j)) — \i4)^ /A. In other cases such as the 
quadratic potential, TZ evolves during inflation due to the presence of isocurvature perturbations, although its change 
is typically small. At the end of inflation, TZ takes almost constant value, 7?. ~ Ci. 
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FIG. 1: The evolution of the adiabatic metric perturbation <l>ad = fc"^"$ad, the 
isocurvature metric perturbation $iso = k'^'^^ieo, and the curvature perturbation 
■^ = k^'^TZ for the case of /3 = 0.09 in the JBD theory with quadratic inflaton 
potential, V{4>) = m^</>^/2. We choose the initial values of scalar fields as cj> — Snipi 
and X = ~l-2»npi, in which case x is nearly zero at the end of inflation. 

One may worry that the above results are obtained by imposing slow- roll conditions, which are only approxim ation s 
t o th e full equations of motion. In order to ans w er s uch suspicions, we numerically solved full equations ( 3.19| )- 
( |3.23 ) along with the background equations (3.1)-(3.4). We adopt the quadratic potential V{(j)) = m^ijP /2, and 
start integrating from about 60 e-folds before the end of inflation. We found that the evolution of field and metric 
perturbations are well described by analytic estimations except around the end of inflation. The evolution of $ for 
the adiabatic and isocurvature mode is plotted in Fig. y, which shows that the contribution of the isocurvature mode 
is small as estimated by Eq. ( 4.15| ). The adiabatic growth of the gravitational potential terminates for rat > 20, 
after which the system enters a reheating stage. During reheating no additional growth of super-Hubble metric 
perturbations occurs in this scenario, unless some interactions between inflaton and other field, a, such as g'^(jP'a'^ /2 
are not introduced. 

The evolution of TZ for j3 = 0.09 (corresponding to ijj — 500) is plotted in Fig. 1, which shows that the change of TZ 
is small during inflation. We have also confirmed that the variation of TZ is negligible in the case of /3 < 0.034. Since 



its growth is sourced by the e''^/'^^'^^ term in Eq. ( 4.1£ ), we find in Fig. 1 that TZ approaches almost constant value 
around the end of infiation. During reheating TZ is conserved except for the short period when H passes through zero. 
Thus one can use analytic expressions for adiabatic curvature perturbations based on slow-roll approximations in 
order to constrain the model parameters of the potential. Note that the number of e-folds, Zk^ after the comoving 
wave-number k leaves the Hubble radius during inflation satisfies 



k 



k 



f 



e-(i-^'/2)^M2zfc)"/2, for 1 « z^ < 60. 



(4.20) 



We can then express the amplitude of curvature perturbation on scale / = 27r/fc as |2£ 

m 



1 



3(1 + w) 
2 



x/2fc^(|CiP) 



27r 



3(1 +w) 



-i^2g/3^.,/2 K^^ f4nzkY 



271 



6n V 



\j n (3 



(4.21) 



which is valid until the second horizon crossing after infiation and also for 1 ^ Zfc < 60. Here w is the ratio of the 
pressure to the energy density. 
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Since the large-angular-scale anisotropy of background radiation due to the Sachs- Wolfe effect is given by 5T /T = 
<I>/3, we can normalize the value of A2n by the COBE normalization [Q. For (3 = 0.034, since this gives the relation, 



^z,^m)-^ --'-^'= ^"v'"pi 



5T _ 1 

= 1.1 X lO"^ 






for n — 1 
for n = 2 



one finds 



m = 1 X 10^^ GeV, n = 1, 
A4 = 2 X 10"", n = 2, 



(4.22) 
(4.23) 



(4.24) 
(4.25) 



which is not much different from the values obtained assuming the Einstein gravity [ plj . Since the behavior of the 
system approaches to that in the Einstein gravity as we increase w, we can conclude that in Brans-Dicke theory, model 
parameters of chaotic inflation should take the same value as in the Einstein gravity. 
Next we consider new inflation with a potential 



V{4>) - ^0 - \c^\ 



for which we flnd, 



V{ip) 



dip = 



^/, 



2A02 



(4.26) 



(4.27) 



;^^ V'iip) " 2A \^02 
We can again express the amplitude of curvature fluctuation as a function of Zk, which is now related with k as, 

(4.28) 



A - „-(l-/3V2)^fc 

kf 



m 



3(1 + u;) 



J^^./i^^2z,f/2 + nHf^^il - e-'^^-/^) 



(4.29) 



with Hf = ^k'^Vo/S. Taking /3 — 0.034 again, it predicts the amplitude of 6T/T to be compared with COBE data as 

ST , , r- Hf , , 

-—(zk ~ 60) ~ 30VA + 0.21^^. (4.30) 

1 rupi 

Since Hf should also satisfy 

(4.31) 



^ < lo-^ 

rripi ^ 



to suppress long- wave gravitational radiation of quantum origin ||62| , we obtain 

A < 1 X 10"", 
from Eq. ( 4.3C|) . Again its amplitude is practically no different from the case of the Einstein gravity. 



(4.32) 



B. Non-minimally coupled scalar field case 

Let us first briefly review the single-fleld inflationary scenario with a non-minimally coupled scalar fleld (^i?0^/2). 
In chaotic inflationary models, Futamase and Maeda |p3] found that the non-minimal coupling is constrained as 
1^1 < 10~^ in the quadratic potential, by the requirement of sufficient amount of infiationn. In the quartic potential. 



This constraint is loosened by considering topological inflation, see 
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such a constraint is absent for negative ^, and as a bonus, the fine tuning problem of the self-coupling A in the 
minimally coupled case can be relaxed by considering large negative values of ^ l|65H 67|. Several authors evaluated 
scalar and tensor perturbations generated during inflation [|68H71[ and preheating | |7^ in this model. Since the system 
is reduced to the single-field case with a modified inflaton potential by a conformal transformation, we cannot expect 
nonadiabatic growth of 7?. on large scales. 

We shall proceed to the case of the non-minimally coupled ip field with Eq. ( ^.10 ) in the presence of inflaton, cj). In 
this theory the evolution of field and metric perturbations was studied in p7| in the Jordan frame. It was found that 
TZ can grow nonadiabatically during inflation on super-horizon scales for negative ^. Here we will show that similar 
results are obtained by the an alysis in the Einstein frame. 

From Eqs. (|4.5| ), ([iq), and ( 1.12 ), we obtain the following explicit solutions: 



Sx = -4(Ci - Cs) 



C^ 



VI - (1 - 60e«2V2 



S(l> = - 



K^V{(I)) 



[Ci{l 



^n'i^') 



■C3(^«V)], 



(4.33) 



* = -Ci-^-C3[6v, 



{^^'^')e4 



(4.34) 



n = Ci-C3 



e-0 - i£.n ^ )f-<p 
e^0 + (l-^K2^2)g^ 



with 



(4.35) 



m^^i^? 



1 - (1 - 60^k2V'2 



1 (V'{d,) 
* 2^2 i^ V{4)) 



(4.36) 



When ^ is negat ive, t he coefficient of -0 in the rhs of Eq. ( 3.27 ) is always positive, which leads to the rapid growth of 
ij} (and x). Eq. (4.33) indicates that long wave 5x fluctuations are amplified with \i}}\ being increased. This is due to 
the fact that the effective mass of 5x becomes negative after horizon crossin g [[l7| , whose property is different from the 



JBD case. In the JBD case, 5x is almost constant during inflation [see Eq. (4.5) with Vi = e '^''x], which restricts the 



nonadiabatic growth of large scale metric perturbations. In contrast, in the present model, <& and TZ exhibit strong 
amplification due to the excitation of lo w mom entum field perturbations unless |'0| is initially very small. 

The second terms in Eqs. (4.34) and (4.35) appear in the presence of non-minimal coupling, whose contributions 
are negligible when j^JK^^s ^ ^ With the increase of \ip\, however, isocurvature perturbations are generated during 
infiation, which can lead to nonadiabatic growth of <I> and TZ. When the second term in Eq. (4.34) grows to of order the 
first term, the adiabatic mode includes the isocurvature mode partially. In thi s cas e one cannot completely decompose 
adiabatic and isocurvature modes in the final results by the expression, Eq. (4.34). 

Let us consider the massive chaotic infiationary scenario with initial conditions, (/)* = imp\ and x* = lO^^TTipj. We 
plot in Fig. 2 th e evol ution of field perturbations for ^ = — 0.02 h i two cases, i.e., (i) solving directly the perturbed 
equations, ( t3.1£ )-( 3.23| ), (ii) using the slow-roll solutions, Eq. ( 4.33| ). In spite of the rapid growth of field perturbations, 
slow-roll analysis agrees reasonably well with full numerical results. The enhancement of Sxk fiuctuations stimulates 
the amplification of S(j)k fiuctuations for mt > 10. Inflationary period ends around mt ~ 17, after which the slow-roll 
results begin to fail. In Fig. 3, the evolution of <& and TZ is depicted. We also plot the first and the second term in 



.(-) 



.(-) 



the rhs of Eq. (4.34), where we denote ^l and $2 : respectively. Although $ is dominated by the ^l term in the 



(^) 



,(- 



,(-) 



is) 



initial stage, ^2 catches up ^\ around mt ~ 7 , after which the ^2 term completely determines the evolution of 
$. We find in Fig. 3 that slow-roll approximations are valid right up until the end of inflation. The growth of metric 
perturbations stops when the system enters a reheating stage. 

If we use the decomposition of Eq. (4.11), the second term gives negligible contribution to the gravitational potential 
around the end of inflation. In the early stage of inflation, however, its contribution is comparable to the first term, 
in which case the isocurvature mode cannot be completely separated from the adiabatic mode. When fluctuations are 
sufficiently amplified, it is inevitable that both adiabatic and isocurvature modes mix each other with the growth of 
fiuctuations, which means that complete decomposition is difficult during the whole infiationary stage. 
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FIG. 2: The evolution of background fields, (j) = "T^/^t-pI and x = x/"^pi) ^.nd large- 
scale field perturbations, S(f> = k'^''^54>/mp\ and 5x = k^ ^x/mpi, in the massive 
chaotic inflationary scenario with a non-minimally coupled x field for ^ = —0.02. 
The initial values of scalar fields are chosen as 0« = Srripi and x* = W~^mpi. The 
slow-roll results of Eq. (i.33) are also plotted, where we denote 5(f)^"' and Sx^"' in the 



figure. We find that slow-roll approximations are valid except in the final stage of 
inflation (mt > 17). The evolution of 5(p^"' and 5x after inflation is not shown. 
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FIG. 3: The evolution of the metric perturbation $ = fc^"<E> and the curvature 
perturbation 7^ = k^^^TZ with same parameters as in Fig. 2. The amplification 
of field perturbations leads to the nonadiabatic growth of $ and TZ. We also plot 
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second terms in Eq. ( f.34 ) , respectively. <1> is mainly sourced by the $2 term, after 



$2 catches up $j^ . The evolution of ^i and <I>2 after inflation is not shown. 
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In Fig. 3 the curvature perturbation, TZ, is nonadiabatically amplified sourced by the second term in Eq. (4.35). 
Whether this occurs or not depen ds upon the strength of the coupHng, ^, and the initial x- When both are small 
and the second term in E q. (4.35) is negligible relative to the Ci term during inflation, we can regard the first and 
second terms in Eq. ( |4.35| ) as adiabatic and isocurvature modes, respectively. In the simulations of Figs. 2 and 3, we 
take X* — 10~^mpi, in which case numerical calculations imply that the conservation of 7?. is violated for ^ < — 0.01. 
When ^ < — 1, strong amplification of TZ is inevitable even for very small values of x far less than ?7ipi. For positive £,, 
conservation of TZ is typically preserved due to an exponential suppression of x during inflation p^ , |7^ , |74| . Regarding 
detailed investigation about the observational constraints of the strength of ^, see [ p7[ whose results are similar to 
those in the analysis of the Einstein frame. 



C. Higher-dimensional theories 



In the higher-dimensional theory with Eq. (2.16), the kinetic term takes a canonical form. In this theory the 
condition, D > 1, gives the constraint, \/2/3 < j3 < -\/2, which is different form the JBD theory with /3 ^ 1. Larger 
values of (3 correspond to the steep exponential potential of the x held, which leads to the rapid evolution toward 
the X direction. In this case inflaton decr eases slowly relative to the x field. Then the expansion of the universe is 
described by the power-law solution, Eq. ( 3.3(]| ). 

For exampl e, in the polynomial inflaton potential, V{(t>) = A2n0^"/(2?^)i classical trajectories of scalar fields are 
given by Eq. (4.2) as 






Differentiating Eq. ( 4.37 ) with respect to t yields 



(4.37) 



Zn n 



2ttP 



mpi 



(4.38) 



This relation indicates that for the values of (j) greater than mpi with P and n being of order unity, |x| is typically 
larger than |0|, in which case x rapidly evolves along the exponential potential. The power-law inflation continues 
until 101 grows comparable to |x|, corresponding to (/>/mpi — n/{\/2nf3). After falls down this value, (j) begins to 
evolve faster than x towar d the local minimum at = in the (j) direction. In this stage the system deviates from 
the power-law ex p a nsion , ( 3.301) . 

From Eqs. (|4.6|) , ( 4.11 ), (^4.12 ), and ( 4.13|) , $ and TZ evolve during inflation as 



$ 



-Ci 



H 



C- 



P^ 



c ^ 






2a/e0 



2(1 + a/) 



(4.39) 



7^ = Ci - Cg 



0" 



13^ + 2e. 



— Ci — C3 



P^ - 2a/e^ 



{l + af){(3^+2e^y 



(4.40) 



where aj = 0^ /{2e^j) is of order unity for ^2/3 < (3 < \f2. Since the condition, (? ^ e^, holds during power law- 
inflation, the curvature perturbation in this stage takes almost a constant value, 7?. ~ Ci — C3 = Ci — C3/(l -t- a/). 



As \4>\ grows relative to |x|, TZ begins to evolve due to the change of e^ in E q. (|4.40| ). This corresponds to the stage 
where deviations from power-law inflation become relevant. When e^ in Eq. (4.40) becomes comparable to /3^/2 (i.e., 
af ~ 1) around the end of inflation, we have TZ ~ Ci — C^af /{I -t- a/) = Ci. After inflation TZ takes this conserved 
value. 



D. The Ft^ theory with a non-minimally coupled x field 



In the f{R) theories, effective potentials do not gen er ally t ake separated forms, Eq. (4.1), as found in Eq. ( 2.21 ). 
Nevertheless we have closed form solutions, Eqs. ( 3.45| )-( p!4£ ), by which the evolution of cosmological perturbations 
can be studied analytically. 
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Let us analy ze the R^ in flationary scenario with a non-minimally coupled x field as one example of the f{R) theory 
[see Eqs. ( ^.22 ) and ( ^.23 )]. Note that when x = the system has an effective potential, 



U = 



'pi 



(327r)2g 



1-e 



-(\/6/3)k(/> 



(4.41) 



The (f> field defined by En . ( ^.19 ) plays the role of an inflaton and leads to inflationary expansion of the Universe in 



the region, cj) > rripi |5|,|75| . In the absence of the non-minimally coupled x field, the resulting spectrum of density 
perturbations after the end of inflation was found in JT^] (using equations for perturbations of the FRW model for the 
Einstein gravity with one- loop quantum corrections derived in [|77|) and then rederived in |78|. Here we study how 
the effect of non-minimal coupling alters the adia batic evolution of cosmological perturbations in the single field case. 
In the presence of non-minimal coupling, Eq. ( p.4q ) is reduced to 



J = 



[1 - (1 - e«'x')e'^^^'^"1(l - C^'x') 
[1 - (1 - C'«2x?)e-(^/3).0,](l _ ^^2^2) ■ 

Then Eqs. (3.43) and ( |3.44 ) are integrated to give 



2^6(1 -e«V)e-(^/')"'^ 
'3K[l-(l-eK2x2)e-(v^/3)«'^] 



Ci — C3 



^^Hx'-xl) 



[1 - (1 - eK2x2)e-(v^/3)«'?i.](l - ^k2^2) 



(4.42) 



(4.43) 



where Ci = 



Sx = - 



-^ "^ ^ ^ l-il-^K^xl)e-^^/'^^^ 

1 - (1 - ^K2;^2)g-(\/6/3)K0 ^ ^ 1 - (1 - ^K2^2)g-(V6/3)K<^. 



4ex 



-K^Qi and C3 = —n'^Q'i- Therefore $ and TZ 



-Ci 



H_ C^ r S,n\x^-X: 

H^ 1 - (1 - C«'x^)e-(^/')"'^* r 1-?«V 



are expressed as 

g^^(x^-x;) 



- e. 



|g(V6/3)K0_(^_^^2^2)| 



(4.44) 



(4.45) 



n = Ci 



C, 



H i^\x^ - Xl) - ^x {e(^/3)''^ - (1 - ^>exl)} 



1 - (1 - ^K2x2)e-(^^/3)«'^* 



£0 + e(\^/3)'^'^ex 



where e^ and e^ are defined by 



1 fu 



2k2 V U 



1 /[/ 



4 / (1 ^ gK2x^)e-(^/3)"'^ \ 
3 I 1 - (1 - ^K2x2)e-(v^/3)'"^ ) ' 



2k2 V U 



S,Kxe 



-(V6/3)ft0 



1 - (1 - ^K2^2)p-(\/6/3)fc0 



(4.46) 



(4.47) 
(4.48) 



In the absence of non- minimal coupling, one finds adiabatic results, $ = —CiH/H^ and TZ = Ci. In two-field inflation 
with a non-minimally coupled x field, the presence of isocuryature perturbations can lead to nonadiabatic growth of $ 
and TZ as found in the second terms in Eqs. (4.45) and ( 4.46| ). Their contributions are negligible when the conditions, 
|^|k;2x^ ^ 1 and e^ <C 1, holds during infiation. The latter condition is similar to the former one when |^| < 1. 



From Eq. (3.27), since x is approximately written as 



A^H 



l-(l-^K2^2)g-(v^/3)K0^' 



(4.49) 



the X field exhibits exponential decrease for positive values of ^. For negative ^, however, x is exponentially amplified 
during inflation, which means that the condition, |C|k^x^ ^ Ij can be violated. The long-wave Sx fluctuation grows 
with th e inc rease of x as found in Eq. (4.44). On the other hand, the growth of S4> begins only when the second term 
in Eq. (4.43|) becomes comparable to the first term. 
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FIG. 4: The evolution of background fields, </!> = (f>/m^i and x = x/nipi, and long- 
wave field perturbations, 50 = k^'^Scfi/m-pi and 5x = k^''^Sx/''npi as a function of 
time, i = mpit/y^96Tvq in the R^ infiationary scenario with a non-minimally coupled 
X field for ^ = —0.025. The initial values of scalar fields are chosen as 0* = l.lrripi 
and X* ~ 10~^mpi. We also plot the slow-roll results, where we denote 5(f)^"' and 
dx in the figure. The evolution of Stp^"' and Sjc"' after infiation is not shown. 
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FIG. 5: The evolution of the metric perturbation $ = fc^"<l> and the curvature 
perturbation ^ = k^'^TZ in the Fl? inflationary scenario with same parameters as in 
Fig. 4. We also plot I- J"' = fc^/^^^J"' and j>f,"' = fc^/^^^"', where $5"' and $<"' denote 



$2 after inflation is not shown. 



the first and the second terms in Eq. (4.45), respectively. The evolution of $j and 



s(=) 



20 



We plot the evolution of 5x, S(j), and the slow-roll results ( 4.43| ) and (4.44) for ^ = —0.025 with initial conditions 



l.lTOpi and X* = 10 ^iripi. In this case inflation ends around t = mp\t/y/96Trq w 130 with e-foldings, N w 63. 
Again the slow-roll analysis is quite reliable except around the end of inflation. 

In Fig. 5 we also depict the evolution of $ and TZ, and the first {— ^i ) and second {— ^2 ) terms in Eq. (4.45). 



In the initial stage of inflation where x 8'iid Sx are not sufficiently amplified, the gravitational potential is dominated 

{ s) ( s) ~ ( s) 

by the ^\ term, in which case $2 ™ay be regarded as the isocurvature mode. However, after i w 70 where <I>2 
catches up $1 , we find in Fig. 5 that <I>2 mainly contributes to the gravitational potential. As is similar to the case 
of the chaotic inflationary scenario with a non-minimally coupled x field, adiabatic and isocurvature modes mix each 
other with the growth of the x fiuctuation. If one defines the isocurvature mode as the one which gives negligible 
contribution to the gravitational potential, it can not be completely separated from the adiabatic mode during the 
whole stage of inflation. 



The conservation of TZ is typically violated when the term proportional to C'3 in Eq. ( 4.46 ) surpasses the one 
proportional to Ci. For the initial values, (/)* = I.lmpi and x* — 10 '^"^pb TZ exhibits nonadiabatic growth for 
^ < — 0.02. Negative large non-minimal coupling such as ^ < — 1 leads to strong amplification of TZ unless x is 
initially very small. Although we do not make detailed analysis here, the basic property is quite similar to the case of 
the subsection B. These results are also expected to hold for other infiationary models with a non-minimally coupled 
X field, since the scalar curvature is proportional to the potential energy of inflaton which slowly decreases during 
inflation. 

V. CONCLUSIONS 

We have studied generation and evolution of adiabatic and isocurvature perturbations during multi-field infiation in 



generalized Einstein theories. Most of the generalized Einstein theories recast to the Lagrangian (2.1) in the Einstein 
f rame by conformal transformations. While behavior of adiabatic perturbations is universal and is given by Eqs. 
(|3.9, 3.11) in the long-wave limit (fc — > 0), isocurvature perturbations behave differently depending on specific gravity 



theories. Making use of slow-roll approximations, we have obtained closed form solutions for all non-decaying field 
and metric perturbations in the long- wave limit. The existence of isocurvature perturbations may lead to significant 
variations of the curvature perturbations TZ and C- 

In this work we considered the following four gravity theories. 

(1) The Jordan-Brans-Dicke theory with a Brans-Dicke field x and inflaton (j). Using the Brans-Dicke parameter 
constrained by observations, the isocurvature mode in the gravitational potential $ is negligible relative to the 
adiabatic mode. Therefore the variation of TZ is typically small in this theory. In particular, for the quartic potential, 
TZ is conserved in the slow-roll analysis. 

(2) A non-minimally coupled scalar field x in the presence of infiaton (p. When the coupling ^ is negative, x and 
its long-wave fiuctuations exhibit exponential increase during inflation, leading to the nonadiabatic amplification of <& 
and TZ due to the existence of isocurvature perturbations. Even in this case we find that slow-roll analysis agrees well 
with full numerical results. When field and metric perturbations are sufficiently amplified, adiabatic and isocurvature 
modes of the gravitational potential mix each other. 

(3) Higher-dimensional Kaluza-Klein theory with dilaton x and inflaton (f>. In this theory the inflationary period 
can be divided into two stages: the first is the power-law infiationary stage where x evolves along the exponential 
potential and the second is the deviation from power-law infiation due to the rapid evolution of ((> around the end of 
inflation. In the former stage TZ is nearly constant but its change occurs at the transition between two stages. 

(4) R^ theory with a non-minimally coupled scalar field x- This system has an additional scalar field playing 
the role of infiaton after conformal transformations. Although this theory has a coupled effective potential which is 
different from the above theories (l)-(3), we have integrated forms of long-wave field and metric perturbations by 
slow-roll analysis, which is found to be quite reliable right up until the end of infiation. Negative non-minimal coupling 
again leads to the nonadiabatic growth of <& and TZ, in which case complete decomposition between adiabatic and 
isocurvature modes is difficult. 

While we analyzed generalized Einstein theories involving two scalar fields, there exist other multi-field infiationary 
scenarios such as hybrid inflation |7|] and models of two interacting scalar fields Q • During preheating after inflation, 
there has been growing interest about the evolution of cosmological perturbations for the simple two-field model with 
potential C/(0, x) — A(/)'*/4 -I- g'^cfPx'^/'^ IO'EII' ^^ ^^ certainly of interest to constrain realistic multi-field infiationary 
models based on particle physics using density perturbations generated during infiation, together with constraints by 
gravitinos p2[ and possible primordial black hole over-production during preheating 
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